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Abstract 

We present a solution of the solar neutrino deficit using three flavors of neutrinos within the 
R-parity non-conserving supersymmetric model. In vacuum, mass and mixing is restricted to the 
v^-Vt sector only, which we choose in consistency with the requirements of the atmospheric neutrino 
anomaly. The v e is massless and unmixed. The flavor changing and flavor diagonal neutral currents 
present in the model and an energy-dependent resonance-induced v e -Va mixing in the sun result 
in the new solution to the solar neutrino problem. The best fit to the solar neutrino rates and 
spectrum (1258-day SK data) requires a mass square difference ~ 10~ 5 eV 2 in vacuum between 
the two lightest neutrinos. This solution cannot accommodate a significant day-night effect for 
solar neutrinos. 
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Neutrino oscillation is the most popular solution of the solar neutrino problem [p], |2| and 
the atmospheric neutrino anomaly ||]. Oscillation in vacuum or in matter, through the MSW 
resonance mechanism, posits that neutrinos have non- vanishing, non-degenerate masses and 
that the basis defined by these eigenstates does not coincide with the flavor basis. 

Supersymmetry (SUSY) with .R-parity non-conservation is an extension of the Standard 
Model (SM) which is consistent with all particle physics experiments and is phenomenologi- 
cally rich || . It carries within it new interactions between leptons and quarks which violate 
baryon (B) and lepton (L) number. In this work we show that the flavor changing neu- 
tral currents (FCNC) and flavor diagonal neutral currents (FDNC) due to the L-violating 
interactions induce mixing amongst neutrinos in matter, the key feature in this alternative 
solution of the solar neutrino discrepancy, even though, in vacuum, the v e state is massless 
and does not mix with the other neutrinos. We also indicate how in this model the parame- 
ters can be chosen to address the atmospheric neutrino anomaly, solving the solar neutrino 
problem at the same time. 

Origins of neutrino oscillation other than mass-mixing, notable among them being non- 
standard interactions of neutrinos with matter, like FCNC, were examined by Wolfenstein 
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It is noteworthy that FCNC and FDNC interactions can drive neutrino oscillations even for 
massless neutrinos without any vacuum mixing through an energy-independent resonance 
effect ||. However, these solutions require large L- violating couplings near their present 
experimental upper bounds ||. This has been examined earlier in connection with the 
solar |6|, [7|, U and atmospheric neutrino data || [H| in the two flavor oscillation framework. 



The new explanation of the solar neutrino deficit that we propose, in contrast, relies on 
an interplay between the ^-interactions for three flavors of neutrinos with matter and their 
masses, keeping vacuum mixing restricted to the — v T sector only. 

Imposing baryon number conservation, we focus on the following L-violating terms in 
the superpotential: 

W = ^ijk Li Lj E c k + X' ijk Li Qj Dl, (1) 

assuming that bilinear terms have been rotated away with appropriate redefinition of su- 
perfields. Here, i,j, and k are generation indices, L and Q are chiral superfields containing 
left-handed lepton and quark doublets and E and D are chiral superfields containing right- 
handed charged-lepton and <i-quark singlets. There are nine A (antisymmetric in (ij)) and 
twenty-seven A' couplings, only a few of which will be relevant for this analysis. 

The interaction of neutrinos with the electrons and <i-quarks in matter induces transitions 
(i) Vi + e — > Uj + e, and (ii) z/j + d — ► uj + d. (i) can proceed via W and Z exchange for i — j, 
as well as via A couplings for all while process (ii) is possible through A' couplings and 
squark exchange. Here we concentrate only on the A-induced contributions. 

The time evolution of the neutrino flavor eigenstates (z/j, i = e,/i,r) is governed by 
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where S± = (m\ ^mfy/AE, T\ = S- cos 2023 V , T 2 = S- sin 29 2 3 V , A\ = y/2Gpn e , A 2 = 
Gf^n/ 'V%, and Rij = \iki^jkin e / '4m 2 '. E is the neutrino energy and 6 2 3 V the vacuum mixing 
angle in the — v T sector, and n e are the neutron and electron number densities in 
matter and rh is the slepton mass. A\ and A 2 in h matter arise from SM charged and neutral 
current interactions, respectively. In vacuum, h matter = and H contains mixing only in the 
Vfi — Vr sector. In h matter , we chooseQ k = 2 in the matter- induced contributions Rij. For anti- 
neutrinos, the time evolution is determined by a similar total hamiltonian H = H° — h matter . 

To obtain the mass eigenstates, first we rotate by U' = (where Uij is the standard 

rotation matrix) and write the effective mass squared matrix, ^ = H — E — A\ — A 2 , in 
the new basis as 

^ 2 ( -Rll c 13 _ 2.R13C23S13C13 + A + S? 3 --R13S23C13 

-^13-S23Cl3 A_ -Rizs^siz I , (3) 
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1 In view of the antisymmetry of \j k in in order to generate the mixing of the v e with the other 

neutrinos we have to choose k = 2 or 3. For the latter choice, mixings due to ft interactions are very small; 
for example, A131A231 is highly constrained from fj, — > 3e decay 0. 
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and Cij = cos 6 ij and = sin 6V,-. Furthermore, 



tan 2^3 = 2X2/(2^ + ^33); tan20 13 = 2i2 13 c 23 /£)i; £>i = A+ - i? n . (5) 

Note that #23 ~ #23,, whilef^] 6 , 13 ~ except near a possible resonance, when Di = 0. We 
show below that this resonance condition cannot be achieved in the sun. Consequently, to a 
good approximation, the third state in this basis decouples in eq. The upper left 2x2 
block is readily diagonalised, resulting in three effective masses rrtj as: 



m\/ (2E) = c\ 2 [Rnc\ z - Ri 3 c 23 sin 26 13 + A+s? 3 ) + i?i 3 s 2 3Ci3 sin 26 12 + As 



m 2 /(2E) = s l2 [Rnc 13 - Ri 3 c 23 sin 26»i 3 + A + s 13 J - i?i 3 .S23Ci3 sin 26 12 + A_c 12 

ml/{2E) = R u sj 3 + R 13 c 23 sin 20 13 + A + c? 3 , (6) 
where 

t ail 29 12 = ~ 2R ^ C ^ . D 2 = A_ - R llC j 3 + R 13 c 23 sin 2 13 - A + sj 3 . (7) 
u 2 

A resonant enhancement of Q\ 2 occurs when D 2 = 0. 

The neutrino flavor eigenstates v a = v e M T are related to the mass eigenstates v% = v\ 23 

by 



E^i, (8) 



where U a i are elements of the unitary mixing matrix 
U = 



^ c 12Cl3 S12C13 S13 

— S 12C23 ~~ c 12 s 23 s 13 C12C23 — S12S23S13 S23C13 (9) 
\ S12S23 ~~ C 12C23 S 13 ~~ C12S23 — S12C23S13 C23C13 



We have chosen real L violating couplings and as such there is no CP violating phase in the 
above mixing matrix. Further, in order to satisfy < 6i 2 < tt/2 in eq. (H) for convenience, 
we take A121A321 < 0. 

As noted above, level crossings and resonance behavior, which are energy dependent 
due to neutrino masses, can occur in two situations, namely, (a) when D\ = 0, and (b) 
when D 2 = 0. Of these, only the latter can be satisfied inside the sun, as we now discuss. 
The sub-GeV and multi-GeV zenith angle dependence of atmospheric neutrinos as well as 
the energy dependence of the up-down asymmetry require Am32 ~ m 3 « 10 -3 eV 2 with 
maximal vacuum mixing in the — v T sector ||. The presence of L violating interactions 
does not alter this significantly (see later). On the other hand, n e at the core of the sun 
is about 1.13 x 10 12 eV 3 . Thus even for E as high as 20 MeV, it is not possible to satisfy 
the (a) resonance condition and hence we consider only the (b) resonance in the subsequent 

2 This follows as A + — R\\ ~ rri^/(2E) is very large with respect to -R13 in the sun. 
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discussion of the solar neutrino data. At resonance, 9 12 = vr/4, while the other mixing angles 
are # 13 ~ and #23 — @23v Recall that away from resonance, #12 ~ and for vacuum 
propagation only 6*23 = 9 2 3 V is non-zero in eq. (||). At first glance, one might think that 
if U13 in vacuum is very small then solar neutrinos will be almost unaffected by the mass 
of v T and analysis with three neutrino flavors may not be essential. However, unlike in the 
SM, where only v e interactions with matter are relevant for neutrino oscillation, in the ft 
super symmetric Model, FCNC and FDNC interactions of all three flavors of neutrinos turn 
out to be important. In fact, one can see from eqs. (|2|) and (0) that R 13 , arising from FCNC 
interactions, appears in tan 2#i2 and plays a pivotal role. 

We now turn to the oscillation of solar neutrinos due to their interaction with matter 
inside the sun. As already discussed, v e in the sun can experience only one of the two 
resonances. S13 in eq. ([]) is very small as noted earlier and we use the survival probability 
of f e valid for a two flavor analysis: 

Pu e ^u e = g + Q - Pjump^j COs2#i2(£l) COS 29 12 (x ) , (10) 

where Xq is the production point inside the sun and X\ the detection point at earthP]. The 
jump probability is Pj ump ~ exp[— 7rj res F/2\, ^ res being the adiabaticity parameter. F = 1 
for the exponential density profile since the vacuum mixing angle is zero and 

rh\ — fh\ m\ f p\ 2 ( n t 



lres AE9 12 ~ E U UvL' 

where k = (2A? al - A| 21 ) /(8mJ) + y/2G F ~ y/2G F and p = \^^\. 

In order to obtain the best-fit values of Ami 2 ~ m 2 and p, we have performed a \ 2 
analysis using the Standard Solar Model (SSM) ]TT[ and the solar neutrino rates from the 
Homestake (CI), Gallex, Sage, and Kamiokande (K) experiments 0. We have also used the 
latest SK rates and spectrum data for 1258 days ||12|| . Taking into account the production 



point distributions of neutrinos from the different reactions (e.g., pp, pep, 7 Be, 8 B etc. 
we have calculated the averaged survival probabilities using eq. (0). Here, we include a 
parameter Xb to take into account a possible deviation of the overall normalization of the 
8 B flux from its SSM value. We set 9 2 3 V = tt/4. The best-fit values of the parameters are 
presented in Table 1 along with Xmim ^ ne goodness of fit (gof), and the calculated rates 
using these values of the parameters. Case (1) is a fit to the total rates. Note that the 
best-fit parameters result in an unusually good fit to the CI rate and the Ga prediction is 
right near the average of the Sage and Gallex data. We have found that the fit improves 
even more if the K-rate is excluded. In case (2) we have fitted the SK spectrum while (3) is 
a fit to the total rates and the SK spectrum^. In Fig. 1 is shown the calculated spectrum for 
SK for the best-fit parameters along with the experimental data. Also shown is the charged 
current spectrum expected at SNO for one sample case, the best-fit values in case (3)Q. 

3 Notice that cos2#i2(xi) = 1, corresponding to 6*12 = in vacuum. 

4 Ai2i and A321 are tightly constrained 0. Besides significant cancellation between these terms is possible 
if they are of same order. 

5 We have dropped a small contribution from the hep process. 

6 We have checked that the fit (3) is essentially unchanged if the SK rate is excluded from the fit. 
7 For these best- fit values, the prediction for the neutral and charged current rates at SNO, normalized to 
the SSM, are 0.56 and 0.44, respectively. 
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Case 


Best-fit Values 


Corresponding Rates 


P 

(io- 24 

eV~ 2 ) 


m 2 

(10~ 5 
eV 2 ) 


X B 


dof 


x l 

(gof) 


CI 
(0.33± 
0.029) 


Gallex 
(0.52± 
0.06) 


Sage 
(0.60± 
0.06) 


K 
(0.54± 
0.07) 


SK 
(0.451± 
0.016) 


1 


0.595 


1.063 


0.845 


2 


1.71 

(42.5) 


0.326 


0.561 


0.561 


0.478 


0.455 


2 


0.009 


0.01 


0.446 


16 


18.76 
(28.1) 


0.582 


0.947 


0.947 


0.446 


0.446 


3 


0.360 


0.980 


0.560 


21 


25.53 
(22.5) 


0.364 


0.558 


0.558 


0.456 


0.447 



Table 1: The best-fit values of the parameters, p 



4m? 



, m|, and Xb from fits to (1) all rates, 



(2) the SK spectrum, and (3) rates and SK spectrum. The rates for the different experiments 
obtained using these best-fit parameters are also shown. 

The best-fit values of Tjl couplings in Table 1 are consistent with the existing constraints. 
For example, in case (3), choosing ~ 100 GeV, we get Ai 2 iA 32 i ~ 0.0144. Ai 2 i is 
constrained from fi —>■ ez/gZ/^ decay (with selectron exchange tree level diagram apart from 



evv)jY{r — > [ivv) 
200 GeV or more, 



the SM W exchange diagram). The bound on A321 is from R = T(r 
which gets a contribution from a selectron exchange diagram. For m- t 
the requirements are easily satisfied J3J. 

Turning now to atmospheric neutrinos ||, for a simple-minded analysis we can consider 
the earth to be a slab of a single density. n e in earth lies in the range (3 — 6)7^ cm -3 . So the 
resonance condition, D2 = 0, cannot be met for atmospheric neutrinos having energy near 
the GeV range. In order to explain the observed zenith angle dependence, we must choose 
Am 32 ~ 10~ 3 eV 2 . This precludes the occurrence of the other resonance, Di = 0. Since 
neither resonance condition can be satisfied, there will be almost no effect on atmospheric 
neutrino oscillation due to the L violating interactions as the associated couplings are very 
small. So one can consider the mixing matrix in eq. (f|) valid for vacuum for which only 623V 
is non-zero. Thus the solution to the atmospheric neutrino anomaly is just the standard two 
neutrino mass-mixing one. 

The neutrino masses and mixing pattern in vacuum required in this solution can naturally 
arise in many models. For example, the trilinear couplings in eq. (0) contribute to the 
neutrino mass matrix at the one-loop level through slepton or squark exchange diagrams 
||13||. In particular, from the A' couplings one obtains: 



m 



loop 



3 m 2 (Ab + /^tan/3 
8n 2 ml 



A i33 A j33- 



(12) 



where Ajj and are soft SUSY-breaking parameters, m^ is the 6-squark mass and tan f3 is the 
ratio of two Higgs vacuum expectation values. The last two generation indices in A' have been 
chosen as 3 for which the loop contributions are enhanced via the b— quark mass. We remark 
that m,ij is very small when i = 1 and/or j = 1 because of the more stringent constraint || 
on A' 133 . Notice that this mass matrix can correspond to almost maximal mixing for and 
if A 233 ~ A 333 , with two neutrino masses very small and one neutrino having significantly 



higher mass 7713 w 2 m 33 op , which can be suitably chosen by taking appropriate values of the 
different parameters in (|T2"|). It should be borne in mind that m 2 depends on the difference 
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of A 233 and A333 and can be several orders less than the mass of the heavier neutrino while 
there will be almost maximal mixing. The remaining neutrino mass is m\ ~ 0. Thus masses 
and vacuum mixings can be as required in the model under consideration. 

This neutrino mixing pattern also satisfies the bound U\^ 2 < 0.04 in vacuum from the 



'13 

0. 



CHOOZ reactor experiment [II]. In fact, in vacuum U13 2 

A comment about the earth regeneration effect for solar neutrinos is pertinent. The v e is 
unmixed with the other neutrinos in vacuum. As n e in earth is about two orders less than 
that near the core of the sun, no resonance condition will be satisfied^. Hence, there will not 
be an earth effect for solar neutrinos. In comparison with the small angle MSW fits [ 15| , the 
somewhat larger best-fit Ami 2 and the zero value of 6*12 in vacuum here, result in a smaller 
day- night effect. 

Though our discussion has been within the framework of i?-parity violating SUSY, there 
are other models |16] where FCNC and FDNC interactions are present. Our results can be 
adapted to these scenarios in a straight-forward manner. 
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Figure 1: The calculated SK solar neutrino spectrum for the best-fit parameters Ami2, p, and Xb 
from (1) rates, (2) SK spectrum, and (3) rates and spectrum. The SK 1258-day data [12] and the 
predicted SNO charged current spectrum for the best-fit (3) are also shown. 
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